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AERODYNAMIC PROPERTIES OF ROUGH 
SPHERICAL BALLOON WIND SENSORS 


SUMMARY 

A first-order theory of the fluctuating lift and drag coefficients associ- 
ated widi the aercH^namic^ly induced motions of rising and falling spherical 
balloon wind sensors is developed. The equations of motion of a sensor are 
perturbed about an equilibrium state in which the buoyancy force balances tbe 
mean vertical drag force. It is shown that to within first-order in perturba- 
tion quantities the aerod 3 mamic lift force is confined to the horizontal, and the 
fluctuating drag force assoc^ted with fluctuations in the drag coefficient acts 
along the vertical. The perturbation equations are transformed with Fourier- 
StielQes integrals and the resulting equations lead to relationships between the 
power spectra of the aerodynamically induced velocity components and the 
spectra of the horizontal lift and drag coefficients. 

The aerod 3 mamically induced motion of the Jimsphere balloon occurs 
predominantly in the horizontal plane, so that ^Lj.ms^^^Drms’ ^Lrms 

and denote the root mean square horizontal lift and vertical drag 

coefficients. The aerodynamically induced motion is an extremely narrow- 
band process with essentially delta function behavior characterized by fre- 
quency cu (radsec"^). The nondimensional frequency co D/w (Strouhal 

number) and variance (o^/w) of the induced zonal or meridional components 

of velocity are functions of the Reynolds number based on the diameter D and 

the mean rise rate w of the sensor. The experimental range of variation of the 

Reynolds number is 1.4x 10®to 6. 6 x 10®. The calculations show that 

Ct /C_ “ 0. 36, where C_ is the mean, or zero order drag coefficient. 
Lrms/ D D ® 

The theory shows that the first-order Fourier components of the 
horizontally induced lift force lead the horizontally induced velocity vector, 
and that the induced vertical drag lags the induced vertical velocity for rising 
balloons and leads for falling balloons. The phase angles of the induced lift 
and drag associated with frequency o> of the Jimsphere balloon are functions 
of the Reynolds number. 



INTRODUCTION 


In the continuing pursuit of hi^ quality wind data derived from meas- 
urements of balloon motions, three factors must be considered: (1) the re- 
sponse of the wind sensor to the environmental wind, (2) the accuracy of the 
tracking system, and (3) the aerodynamically induced motions of the sensor. 
This report is concerned with the tMrd factor. Aerodynamically induced bal- 
loon motions produced by aerodynamic lift forces normal to the velocity vector 
of the air relative to the balloon are the result of the shedding of vorticity. 
Clearly, these induced motions introduce errors in the wind calculations. 

The physical characteristics of the balloon (shape, wei^t, distribution of 
mass, etc. ) and the dynamic and thermodynamic properties of the air are the 
factors that control the induced sensor motion. To filter the induced motions 
from the balloon position measurements in an optimum manner, it is necessary 
that their spectral properties be known. These spectral properties can be 
specified in terms of the spectra of the aerodynamic lift and drag force coeffi- 
cients. Once the power spectra of the aerodynamic coefficients are known, 
it is possible, in principle, to develop filters which can remove the aerody- 
namically induced motion from balloon tracking data. Knowledge of the coeffi- 
cients will also aid in the design of meteorological balloons. 

The purpose of this report is twofold: (1) to develop a linear pertur- 
bation theory of balloon motion, whereby it is possible to calculate spectra 
of the aerodynamic coefficients from self-induced balloon velocity spectra, 
and (2) to apply this theory, to extend our present knowledge of the aerody- 
namic properties of roughened, spherical, balloon wind sensors like the 
Jimsphere. 


The Jimsphere, a roughened balloon two meters in diameter, is widely 
used in the aerospace and meteorological communities to obtain accurate, 
high resolution measurements of the wind profile in the first 18 km of the 
atmosphere. An internal pressure of about 5 mb is maintained to insure con- 
stant volume. The 398 conical roughness elements, each approximately 7. 5 
cm wide at the base and 7. 5 cm high, serve to control vortex shedding in the 
supercritical Reynolds number region below 11 km; i. e. , they decrease the 
spectral bandwidth of aerodynamically induced motions. 
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BASIC EQUATIONS 


The basic equations of motion which govern the motion of a balloon or 
falling sphere wind sensor are given by the vector equation 

-Vl» 

“a* 

where m, m , A, C,. , and denote the mass, apparent mass, cross- 
a L D 

sectional area, the lift coefficient vector, and drag coefficient of the sensor, 
respectively. The symbol m^ is the mass of the air displaced by the sensor, 

g is the acceleration of gravity, p is the density of the atmosphere, and t is 

the time. The quantities V and denote the velocity vectors of the sensor 

and the environment, and are given by 

V = ui+vj+wk (2) 

V^=ui+vi+wk , (3) 

e e e •" e 

where i, j, and k denote unit vectors in the zonal (x) , meridional (y) , and 

vertical (z) directions, and u, v, w and u , v , w denote the corresponding 

0 0 0 

velocity components of the sensor and environment. 

The first and second terms on the right-hand side of (1) represent 
the drag and lift forces. The drag force is parallel and opposite in direction 

to the wind vector relative to the balloon (V - V). The lift force is perpen- 

e 

dicular to the relative wind vector. The third term represents the apparent 
mass effect [ 1] , and the fourth term the buoyant force. If the mass of the 
air displaced by the sensor is greater than the mass of the sensor (m<m^) , 

the balloon will experience a positive buoyant force. Conversely, if m is 

greater than m , the sensor will experience a negative buoyant force. If 
o 

m equals m^, the buoyant force vanishes and the balloon is said to be neutral. 
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FIRST-ORDER PERTURBATION EQUATIONS 


The environmental wind is assumed to be composed of two parts: (1) 
a constant basic state horizontal flow , and ( 2) a superimposed velocity per- 
turbation., so that 


u (t) = u + 

u'(t) 

e e 

e 

V (t) = V + 

v’(t) 

e e 

e 

CD 

11 

CD - 



The overbars and primes denote the basic state and fluctuating parts of the 
wind. Similarly, the velocity componehts of the sensor are represented by 

u(t) == u -I- u'(t) I 

v(t) ^ V + v'(t) > . ( 5 ) 

w(t) = w -F w*(t) I 

Unlike the assumed environment, the sensor has a mean vertical velocity equal 
'to w. One assumes that the perturbation quantities in equations (4) and (5) are 
infinitesimal or sufficiently small so that second- and hi^er-order terms in 
perturbation quantities can be neglected with respect to first-order terms. 
Furthermore, it is assumed that the horizontal mean motion of the sensor is 
in equilibrium with the mean flow of the environment; thus , 

u = u 

e 

V = V 

e 



The Lift Forces 

The aerodynamic lift forces which act in a plane perpendicular to the 
relative wind vector are a result of vortex shedding and instability of the wake. 
The magnitude and direction of these forces vary in time and are responsible 
for the aerodynamic ally induced motions exhibited by wind sensors like the 
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ROSE [2,3] and tiie Jimsphere [4] balloons. The lift coefficient vector is 
given by 


= 


Lx 




+ k 
Lz 


( 7 ) 


where C, 


and Cj are the lift coefficients of the x, y, 

j—iZ 


and z compo- 


'Lx’ Ly’ 

nents of the lift force vector. The coefficients are functions of time and 
will be treated as perturbation quantities. As mentioned previously, the lift 
vector acts in a plane perpendicular to the relative wind vector, and the orien- 
tation of this plane is defined by the following unit vectors which are parallel 
to the lift force plane: 


Tx(V - V) 
e 

lTx(v - v)i 

e 


7x(V -V) 
ljx(V^-V)l 

W 


( 8 ) 


(9) 


Substitution of the perturbation expansion equations ( 4) and ( 5) into equations 
(8) and (9) and utilization of equation (6) yields, to within first-order, 



( 10 ) 


( 11 ) 


These results show that the plane of the lift force is tilted slightly from the 
horizontal plane. The lift coefficient vector can be written as 

The components of this vector can be related to the components as given by 
equation (7) by combining equations (10), (11), and (12), to yield: 
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( I f . .1 

V - V \ /u - u 

W / -^ \ w > 


( 13 ) 


Comparison of equation (13) with equation (7) shows that 

and that is therefore a second-order perturbation quantity, so that the 

Lz 

first-order lift coefficient vector is given by 


i + j 
L Lx Ly 


(14) 


T his result shows the first-order lift coefficient vector to be in the horizontal 
plane. 


Substitution of the perturbation expansion equations (4) and (5) into 
the lift force term in equation (i) and utilization of equations (6) and (14) yield 
the first-order lift force 




1 —9 

F,. = — p A w^ { Ct i + 3 

L 2 V Lx Ly •’ 


(15) 


The Drag Force 

The drag coefficient can be represented in terms of a mean value Cj^ 
and a superimposed perturbation C^: 




( 16 ) 


Substitution of equations (4) , (5) , and (16) into the drag force term in 
equation (1) yields, to witiiin first-order in perturbation quantities, 

= - 4" P AC^Iwlw k + 4" p AC^Iwl Ru' - u‘) 1 + (v' - v') j 


+ 2(w* - w’) k 
e 




— p AC^Iwlw k 


(17) 
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The first term on the right-hand side of equatioii (17) is the mean drag force. 
The second term is the contribution to the drag force resulting from perturba- 
tions in the relative wind vector, and the third term is the contribution result- 
ing from fluctuations in the drag coefficient. 


The Perturbation Equations 

Substitution of equations (4), (5), (15), and (17) into equation (1) 
yields the linear perturbation equations that govern the balloon motion: 

“ T p A 

(18) 

dt 2 D e 2 Ly a dt e 


(19) 


= pAC^Iwl (w' - w*) - pA wlwl C' + m ^ (w' - w*) 

Qt 13 0 2 J-) 3. Clt/ 0 


— pACj^lwlw - (m - m^) g 


( 20 ) 


In deriving equations (18) , (19) , and (20) , it was assumed that the environ^ 

mental wind field has negligible horizontal and temporal variations; i. e. , V* 

0 

is a function of altitude z only, so that the derivatives of the environmental 
wind components with respect to t in equations (18) , (19) , and (20) are the 
changes in the environmental wind as would be seen by an observer attached 
to the balloon ( see Reference 5 for more details) . 

According to the theory of linear perturbation expansions , terms of 
like order balance , so that 


— p A |w |w 


= (m^ - m) g 


( 21 ) 


which permits equation (20) to be expressed as 

dw ^ — — 1 — — * d 

m -7— = pAC |w| (w' - w*) - — pAwlwlC* + m — (w* - w’). 
dt D © 2 D a at © 

( 22 ) 
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The basic equations (18) , (19) , and (22) will be used in the analysis that 
follows. In these equations, the lift forces are confined to the horizontal 
plane , and the fluctuating portion of the drag force resulting from fluctuations 
in the drag coefficient is confined to the vertical direction. 

FOURIER-STIELTJES REPRESENTATION AND 
rms AERODYNAMIC COEFFICIENTS 


The time -dependent quantities in equations (18) , (19) , and (22) can be 
represented with Fourier-Stieltjes integrals; i. e. , 

^ (t) = J°° e^^^dZ (co) , (23) 

— OO ^ 


where |(t) is a time-dependent quantity, and dZ^(co) is the complex Fourier 
amplitude of ^(t) at frequency co and is given by 


dZ (co) 


lim 1 r o^.,. 
O -T 

O 


-icot 

e 



dt 


(24) 


The equations that govern the Fourier amplitudes of the time -dependent 
quantities in (18) , (19) , and (22) are obtained by substituting their Fourier- 
Stieltjes representations and noting that the functions exp(io)t) constitute a 
complete orthogonal function set, so that 


dZ (co) 
u 


dZ (co) 

V 


dZ (co) 
w 


1 + ijucoT 
1 + icoT 


1 + ipcoT 
1 + io)T 


2 + ijUcoT 
2 + icoT 


dZ (oj) + 
u 

e 

dZ (co) + 

V 

e 

dZ (co) - 
w 

e 


|w| 

dZ (co) 


1 + icoT 

Iwl 

dZ^(co) 


1 + io)T 

w 

dZ (co) 
z 

oil 
0 1 

2 + icjT 


(25) 


(26) 


(27) 
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' 1 


where, dZ and dZ are the Fourier amplitudes of u’, v', w* and 

/ U,v,w U^, Vg, Wg 

u' , V* , w' ; dZ are the Fourier amplitudes of and C* . The 

e e e x.y.z Lx.y D 

quantities T and M are given by 
m + m 


T = 


a 


iwi 


|m -ml 
o 


g 


(28) 


m 


M = 


m + m 


(29) 


The power spectra of the velocity components of the balloon motion can be ob- 
tained by multiplying each of the equations (25) , (26) , and (27) by its com- 
plex conjugate and performing an ensemble average over all realizations; 
thus 


, , , 1 + (mwT)2 ,^^2 

0 (oj) = 0 (io) + — — 

1 + (o;T)2 % 


0^(CO) 


i + (o)T)' 


Re 




0 (w) - pcoT Im 

'^UgX 




1 + (wT)' 


(30) 


(co) = 


1 + (/nojT)' 


1 + (o;T)^ e 


0 (cu) + — 


w^ 


0 (o)) 

_z 


D 


+ 2-^ 


1 Re 

0 (W)1 

ill 

. V J 


1 + (o)T) 


1 + (0)T)' 


(31) 


9 


I 



<t> (w) 

w 


4 + (HioT)^ 
4 + (coT) 


e 



^^(co) 

4 + (wT)^ 



where the 4>'s with a single subscript are autospectra and the 4>*s with the 
double subscripts are cross-spectra. Thus, for example. 


0 (o))da) = /dZ (oj) dZ* (co)\ 

e ^ e e / 

0 (co)dw = < dZ (co) dZ* (oj)S 

u X \ u X V 

e ' e / 


(33) 

(34) 


where ♦ denotes complex conjugation, and < > denotes the ensemble average 
operator. The notations Re(0) and lm(0) represent the real and imaginary- 
parts of 0. 

The first terms on the right-hand sides of equations (30) , (31) , and 

( 32) are the contributions to the balloon motion spectra by the environmental 

wind fluctuations. These terms have been discussed in detail in Reference 5. 

The last terms are the contributions resulting from correlations between the 

aerodjmamic coefficients and the environmental wind fluctuations. It is 

assumed that these cross-spectral terms will vanish. The second terms are 

the contributions to the balloon motion spectra resulting from fluctuations in 

the aerodynamic coefficients. These terms correspond to the aerodynamically 

induced balloon motions. Representing these terms by 0 (co) , one can 

u,v,w 


c ^ r 1 

1 + (coT) 2 0 (w) , (35) 


0 (w) 
y 

II 

^1 o 

1 + (coT)^ 


y 

(36) 

0^(w) 


[4 + (o)T)^ 


• 

(37) 
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The lift force is statistically isotropic in the horizontal plane; therefore, 

4> ico) = <p (co) . (38) 

X y 


Thus , only one component of the lift force vector in the horizontal plane need 
be considered. The root mean square (rms) lift and drag coefficients can be 
obtained by integrating equations ( 35) and ( 37) over the domain — <»< ox «> 
with the result 


C 


2 

Lrms 



r 21 

f 1 + (cuT) (co)da) 


C 


2 

Drms 






$ (a))dco 
w 


, (39) 


(40) 


where cr and a are the standard deviations of the zonal and vertical 
u w 

self-induced velocity components. The standard deviation of the meridional 
component is equal to cr^. The quantities #^(co) and #^(o)) are the 

corresponding normalized spectra referenced to the half -interval o < w < °° , 


^ (cu) 
u 


# (co) 
w 


ip^ico) 

cr 2 
u 


^ (co) 
w 


w 


(41) 


Equations (39) and (40) can be written in terms of integrals over the half- 
interval o<cjo <°° because the spectra of the aerodynamic coefficients are 
even functions of co. These equations facilitate the calculation of the rms 
aerodynamic coefficients from relatively simple measurements of the aero- 
dynamically induced sensor motion. 
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AERODYNAMIC COEFFICIENTS FOR 
NARROW BAND PROCESSES 


The self-induced balloon motions of rough spherical sensors like the 

Jimsphere are high frequency, narrow band processes. Figure 1 contains 

spectra of the zonal, meridional, and vertical components of velocity of a 

Jimsphere wind sensor for the altitude interval 4. 4 < z < 5. 4 km calculated 

from FPS-16 radar data taken between 1400 to 1500 u. t. on December 23, 1964. 

The zonal and meridional spectra have peaks at approximately 0. 21 Hz. These 

peaks are attributed to the self-induced balloon motion which results from the 

fluctuating, aerodynamic lift force. The absence of a pronounced spectral 

peak in the vertical velocity spectrum at, or close to, 0. 21 Hz precludes the 

existence of a significant vertical component of self-induced balloon motion 

and implies that « C,, . This means the self-induced balloon 

Drms Lrms 

motion is confined primarily to the horizontal. If the spectral peaks in Figure 
1 are associated with a self-induced spiral motion as observed by Scoggins 
[ 4] , then the zonal and meridional components should be out of phase by 90 
degrees. To verify that this is indeed the case, the Jimsphere velocity data 
were filtered with a Martin-Graham [6] band pass filter, the response function 
of which is shown in Figure 2. The example of the filtered data given in Figure 
3 clearly shows that the zonal and meridional balloon velocities are out of 
phase by 90 degrees. The filtered vertical velocity data, also in Figure 3, 
show that the vertical self-induced motion, if it exists, is very small. 

A particularly useful way to represent a narrow band process is with 
the function 


^u^^^ 1 - F(- o) /e) 


-(oj - oj ) 2 

o / 


€ (27T) 


V2 


, o<oj<°o , (42) 


where 


F(-w A) = 1 - r 
o 


1 

e . 


(27T) 


V2 


(43) 
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Filter Response 



Frequency, Hz 


Figure 1. Power spectra of the zonal, meridional, and vertical components 
of velocity of the Jimsphere balloon observed at 1300 u. t* on December 
23, 1964 for the height interval 4. 4 < z < 5. 4 km. 



Frequency, Hz 


Figure 2. The response function 
used to filter the Jimsphere velocity 
data to isolate the aerodynamically 
induced balloon motion. 


The maximum value of (co) occurs 

u 

at cu = CO , and e is a measure of the 
o 

width of the spectrum. If e Oj then 

(co)^ 6 (co - CO ) , (44) 

u o 


where 6 (o) - oj ) is the Dirac delta 
o 

function [7], 


Substitution of equation (42) into 
equation (39) yields 

t 2 [2 e coJ( -co^/e) - 

1 - F(-cu /e) 

° (45) 
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I -1 - 1 - 1 1 

330 340 350 360 370 

Elapsed Time From Release, sec 


Figure 3. Example of isolated aerodynamically induced Jimsphere balloon 
velocity components. (The zonal and meridional velocity components are 
out of phase by 90 degrees. The vertical component of the induced 
velocity vector is extremely small.) 


where 


1 

fU) = St e (46) 

(27r) 

and f'(-co /e) denotes the derivative of f(£) with respect to £ and 
o 

evaluation at £ = -co /e. If the widths of the balloon motion spectra are suffi- 

o 

ciently narrow, such that 3 say, the quantities fi-co^/e) and 

/e) approach zero, so that the second term on the right-hand side of 
o 

equation (45) can be neglected in comparison to the first term. Thus, the 
approximate expression for the rms lift coefficient is 
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cj 

Lrms 


(47) 





1 + 


T* (co ^ + 
o 




Furthermore, if the balloon motion spectrum is sufficiently narrow, 

such that €«cjo , then 
o 


C 


2 

Lrms 



1 + (w T) 
o 


(48) 


This corresponds to the case in which ^^(co) is a Dirac delta function as 
given by equation ( 44) . 


Even though the magnitude of the vertical self-induced balloon motion 
is extremely small , an estimate of an upper bound on can be obtained with 

the data in Figure 3. The result that corresponds to equation (47) for the rms 
drag coefficient is 


C 


2 

Drms 


a 

- 2 _VJ_ 

D ^2 


4 + T^ (co ^ + e^) 
o 


(49) 


The ratio between the rms drag and lift coefficients is 


Drms 

^Lrms 


a 

w 


a 

u 


4 + T^ (co ^ + 6^) 
o 

1 + T^ (o) ^ + e^) 
o 




(50) 


The ratio is approximately equal to 0. 1, at most, according to the 

data in Figure 3. In addition, an upper bound value of the coefficient of 

cr /a in equation (50) is 2: therefore, ^ 0. 2 .In all pro- 

w/ u ^ Drms Lrms ^ 

babilily, the radar error contained in the w* time history in Figure 3 accounts 

for at least 50 percent of the signal; thus, a more reasonable estimate is 

s 0. 1 , which seems to verify the earlier conjecture that 

Drms Lrms •* 

. The radar error is small in comparison to the signals 

Drms Lrms 

associated with horizontal components of the aerodynamically induced motion. 
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JIMSPHERE LIFT COEFFICIENT 


To determine the rms lift coefficient of the Jimsphere wind sensor, 
three dual radar tracks of the Jimsphere, tracked with the FPS-16 radar, 
were selected for analysis. The measurements were made between 1300 to 
1500 u. t. on December 23, 1964, December 18, 1967, and March 19, 1968, 
at the Kennedy Space Center, Florida. These tests will be referred to as 
flights 1, 2, and 3, respectively. The measurements extended over the al- 
titude interval 0. 5 to 16 km. The radar provided range and elevation and 
azimuth angle data on the location of the balloon at 0. 1-sec intervals. These 
data were converted to Cartesian coordinates and corrected for the curvature 
of the earth. The Cartesian coordinates were differenced over 0. 2-sec in- 
tervals to yield estimates of the zonal, meridional, and vertical velocities 
of the balloon at 0. 1-sec intervals. The resulting velocity time histories 
were used to calculate the rms lift coefficient of the Jimsphere balloon. 


As a constant diameter spherical wind sensor ascends through the 
atmosphere, the Reynolds number 


Re 


lw Id 

p 


(51) 


decreases because the coefficient of kinematic viscosity p usually increases 
with height while the rise rate |w I usually has a slow decrease with height. 
The quantity D is the diameter of the sensor. Since the aerodynamic proper- 
ties of spheres depend on the Reynolds number, one should expect a depend- 
ence of the lift coefficient on altitude. Accordingly, the balloon velocity time 
histories were divided into 1-km segments to calculate the rms lift coefficients 
of the Jimsphere. 

Power spectra of the zonal, meridional, and vertical components of 
velocity were calculated for each time history segment. This resulted in 15 
spectra for each component of velocity for each balloon track, or a total of 
270 spectra. The vertical velocity spectra looked like the example shown in 
Figure 1 , and thus did not reveal any obvious vertical aerodynamically induced 
motion. It was assumed that the zonal and meridional self-induced velocity 
spectra are of the form given by equation (42). This representation proved 
to be reasonable. Rogers and Camnitz [ 8] analyzed line-of-sight velocity 
spectra of the 2-meter ROSE balloon calculated with Doppler radar measure- 
ments and found that equation ( 42) is a valid representation of the aerody- 
namically induced balloon motion. In the case of the Jimsphere balloon, 

F(-co /e) is approximately 0; thus, one may write 
o 
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-(a; - CO ) 
e o 



( 52 ) 


To develop a model of the aerodynamically induced velocity spectra, it was 

assumed that co $ (oj) is a universal function of co/co . The nondimensional 
o u o 

frequency or Strouhal number^ of the sensor 


S 


CO D 
o 


iw I 


(53) 


and the nondimensional variance 
Reynolds number: 



are universal functions of the 


S = (Re) 


The nondimensional power spectrum is given by 


CO $ (co) 
o u 


1 

(2ir) a 


-( 1 - oo/co ) 
o 
e 



where 


(54) 


(55) 


a 



(56) 


The quantity a is a universal constant. The numerical value of o: was 
obtained by determining that value of co/co^, co'‘'/co^, say, at which co^#^(a)) 

is equal to one-tenth of the peak value; that is, (co*) = 0. 1 * 

0. i [ ( 27 T) ^2 Substitution of into equation (55) yields the result 


1. Usually the Strouhal number refers to the nondimensional frequency 

of the dominant Fourier component in the wake. Here, however, it is 
used to denote the nondimensional frequency of the sensor. 
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(57) 


1 - 

o 

(2 In 1 0)^/2 


The quantity a was calculated for each horizontal spectrum with equation 

(57) . The various values of oi showed no dependence on height and thus 

Reynolds number. The expected value of a = 0. 00906. Figure 4 is a plot of 

the nondimensional spectrum co # (oj) as a function of oj/oj for d = 0. 00906. 

o u o 

2 

The peak frequency and the area cr^ under the aerodynamically 

induced spectral peak were used to calculate the nondimensional functions 

F, (Re) and F„ (Re). Figures 5 and 6 are plots of co /‘Itc and cr ^ as 
12 0 / u 

functions of altitude. The plotted data are the expected values of co ^ j 2 tt and 
(T^ for each 1-km segment of each balloon flight obtained by averaging the 

corresponding zonal and meridional values. A smooth line was drawn through 

each set of data. The smoothed values of co / 27 t and a ^ , the standard pro- 
_ 0 / u 

file of w as a function of altitude calculated by DeMandel and Krivo [ 9] , and 

the standard atmosphere profile of v [10] were used to calculate the non- 

dimensional functions F (Re) and F (Re). These functions correspond to the 

Strouhal number and the nondimensional zonal or meridional velocity variance 
(cTu/w)^ (Figs. 7 and 8). Figure 7 shows that the Strouhal number has a 

typical value equal to 0. 55 with a maximum value approximately equal to 0. 60 
at Re = 2. 5 X 10®. Figure 8 shows that the nondimensional variance is a 
monotonically increasing function of the Reynolds number in the interval 
1. 4 X 10® < Re < 6. 7 X 10® with values ranging between 0. 03 and 0. 092. 

The final parameters needed to calculate the rms lift coefficient are 
the nondimensional time constant 


T |w I 

D 


(58) 


and the mean drag coefficient of the Jimsphere. Fichtl [ 5] calculated 

T as afunctionof altitude for the Jimsphere. Substitution of T into equation 
( 58) with the mean rise rate profile of DeMandel and Krivo [ 9] and plotting 
the resulting values of t as a function of Reynolds number yields the curve 
in Figure 9. Figure 10 contains a plot of as a function of Reynolds num- 
ber. 

In view of the extremely small values of d with respect to unity, one 
may treat the power spectra of the zonal and meridional velocity spectra of the 
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Figure 5. The peak frequency 
0 )^y 2 ir of the aerodynamically 

induced motion of the Jimsphere 
as a function of altitude. 


Figure 4. The nondimensional spectrum 

CD # (cd) as a function of cd/cd 
0 u 0 

for Qi = 0. 00906. 






Altitude, km 




0.085 






Figure 8. The nondimensional zonal 
or meridional velocity variance 

as a function of the 


Figure 9. The nondimensional time 
constant t of the Jimsphere as a 
function of Reynolds number Re. 


Reynolds number Re, 



Figure 10. The mean drag coefficient and the rms lift coefficient 


C, as functions of the Reynolds number Re. 
Lrms 





Jimsphere aerodynamically induced motion as Dirac delta functions. This 
means that equation ( 48) is applicable to the problem and can be written in 
terms of S and t as 


= C ' ^ 
Lrms D ^ ^ 


[l + (St)^] 


(59) 


Substitution of S, (cr^y/w)^, r, and as given in Figures 7, 8, 9, and 10, 
into equation ( 59) yields the rms lift coefficient as a function of Reynolds 


number (Fig. 10). The results in Figure 10 show that the 
curves tend to have somewhat similar shapes and a typical value of the ratio 


and 

Lrms D 


Lrms 


= 36. 


PHASE ANGLES 


To determine the phase angles of the fluctuating lift and drag forces 
relative to the aerodynamically induced balloon motion, one writes the second 
terms on the right-hand side of equations (25) , (26) , and (27) as 


dZ = — e 


dZ 


D 


dZ (oj) 

X 


Iwl ^^1 
= — e 


[l + (o)T) 2] 
dZ (oj) 

y 


[l + (coT)2] 


% 


dZ = 
^1 


1^ 




dZ (<x>) 
z 


4+ (t^T)2] 


V 2 


where 


= tan ^ (-wT) 


7T 


^2 = i" r 


W 


IW I 


+ tan 


' (-‘f) 


(60) 


( 61 ) 


22 



RAD 


The quantities dZ are the Fourier amplitudes of aerodynamically 

U^ ) ) w ^ 

induced components of velocity of the balloon. The quantity is the phase 

angle of the zonal and meridional components of the induced lift force relative 
to the corresponding components of the induced velocity and x« is the phase 


angle of the vertical component of the induced fluctuating drag force relative to 
the induced vertical velocity. These results imply that the induced aerody- 
namic lift force leads the corresponding induced horizontal component of veloc- 
ity. For a rising balloon (w>0) , the induced vertical drag force lags the in- 
duced vertical velocity and vice versa for the falling balloon (w <0) . Figure 11 
contains plots of and x« as functions of cuT for the rising balloon case. 

The phase angles for the falling balloon can be obtained from Fig. 11 by sub- 
stracting it radians from x^. 



Figure 11. The phase angles x^ 

and x„ of the Fourier compon- 

ents of the first-order aerody- 
namically induced horizontal 
lift and vertical drag forces 
as functions of caT fur the 
rising balloon (w>0). For 
the falling balloon (w < 0) , 

7T radians must be sub- 
tracted from x„- 


The Jimsphere induced motion 
is characterized by essentially one 
Fourier component located at co = . 

Thus, one writes equation (61) as 
X^ = tan" -St) 

-1 

X 2 = + tan 

(62) 

The quantities S and t are functions 
of the Reynolds number, so that x^^ 

and Xq ^Iso functions of the 

Reynolds number. Figure 12 gives 
X. and Xo functions of the Reyn- 

olds number for the Jimsphere wind 
sensor based on the data in Figures 
7 and 9. 
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CONCLUSIONS 



Figure 12. The phase angles X^(St) 

and Xo (St) of the Fourier compon- 

ents associated with the first-order 
aerodynamically induced horizontal 
lift and vertical drag forces as 
functions of the Reynolds num- 
ber for the Jimsphere wind 
sensor. 


The aerodynamic properties of 
rough spherical wind sensors have been 
established by using linear perturbation 
expansions of the equations of motion 
of a sensor and FPS-16 radar tracking 
data of the Jimsphere balloon. Integral 
relationships between the first-order 
rms drag and horizontal lift coefficients 
were obtained from the first-order per- 
turbation equations. These relation- 
ships show that the rms lift and drag 
coefficients can be obtained from 
measurements of the mean drag coef- 
ficient, the mean vertical velocity, and 
the time constant of the sensor and 
estimates of the component velocity 
spectra of the aerodynamically induced 
balloon motion. These relationships 
are reduced to rather simple equations 
in the case of sufficiently narrow band 
processes. The mathematical analy- 
sis can be applied to both rou^ and 
smooth sensors. 


Analyses of the FPS-16 radar tracking data of the Jimsphere balloon 

showed that the aerodynamically induced motions are extremely narrow band 

processes characterized by velocity power spectra that have delta function 

behavior, for all practical purposes. Since aerod 3 mamically induced motion 

occurs predominantly in the horizontal plane , it is concluded that the first- 

order rms drag coefficient is at most one-tenth of the horizontal rms lift 

coefficient. These results are supported by the earlier results of Scoggins 

[4]. It was found that the nondimensional induced zonal or meridional velocity 

variance, Strouhal number, and the horizontal rms lift coefficient are functions 

of the Reynolds number. They should also be functions of the mass ratio 

m/m , where m is the mass of the sensor and m is the mass of air dis- 
' o o 

placed by the sensor. To determine how the Strouhal number and nondimen- 
sional velocity variance depend on would be rather difficult because 

many sensor configurations would be required to cover a sufficiently wide 
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range of values of m/m^. The horizontal rms lift coefficient is approximately 
equal to 0. 36 Cj^. 

The precise nature of the mechanism which produces the aerodynami- 
cally induced motion of rough sensors , or smooth sensors , is rather compli- 
cated. It is obvious that the wake is unstable and that vorticity is being shed 
from the balloon resulting in periodic x- and y-directed aerodynamic lift forces. 
However, the nature of the instability which results in the oscillating wake is 
not clear. 

George C. Marshall Space Flight Center 

National Aeronautics and Space Administration 

Marshall Space Flight Center, Alabama 35812 
December 28, 1970 126-61-10-00 
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